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In the context of confining gauge theories we study the flux tube generated by a pair of static 
sources belonging to higher rank representations of the gauge group. Using a simple geometric 
approach based on minimal surfaces describing the world-sheet of the underlying k-string, we 
infer logarithmic broadening of the flux tube as a function of source separation. It turns out that 
the coefficient of the logarithm does not depend on the specific representation of the source nor on 
its N-ality and is universal. Numerical results on a Z4 gauge theory in three dimensions with static 
sources carrying two units of elementary flux compare very favourably with these predictions. 
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It is widely believed that a confining gauge theory near the continuum limit is in the rough 
phase. This means, as explained long ago by Luscher, Miinster and Weisz [|l|], that the colour flux 
tube broadens logarithmically with source separation. This phenomenon can be seen as a conse- 
quence of the underlying string description of the infrared properties of the flux tube. The detailed 
form of the string action is not know, however its IR limit is expected to be a two-dimensional 
massless Gaussian model (free string) made with d — 2 bosonic fields (pj^ ((§1,(^2) describing the 
transverse fluctuations of the world-sheet having a Wilson loop W (C) as boundary. These quantum 
fluctuations give rise to two universal effects: one is a contribution to the vacuum expectation value 
of the Wilson loop (W{C)) which depends only on the shape of C. It leads in particular to the 
Luscher term in the static potential. These kinds of shape corrections have been observed in many 
gauge systems in three and four dimensions for various gauge groups. 

The other universal effect is just the logarithmic growth of the mean square width of the colour 
flux tube. Observing this broadening is very challenging from a computational point of view. As 
a matter of fact, an uncontroversial observation of it has been made only in 3D Z2 gauge model, 
thanks to the efficiency of the Monte Carlo algorithms for its dual, the Ising model [p. 

The question naturally arises whether such a broadening is present also in the k-strings i.e. 
the effective strings describing the confining force between static sources carrying the quantum 
numbers of k copies of the fundamental representation. These strings can be considered as bound 
states of k fundamental strings. It is not clear whether in the IR limit the only relevant degrees of 
freedom are the transverse displacements of a single string or whether we have to add some new 
field describing the possible splitting of the k-string into its constituent strings. Therefore we are 
not entitled to apply blindly the IR gaussian description of the fundamental string to the k-string 
issue. 



1. K-strings and minimal surfaces 

To inspect the width of the flux tube between two static sources in the fundamental represen- 
tation one starts from the quantity 

p ^ {Wf{C)Wf{c))-{Wf{C)){Wf{c)) 

where Wf{C) and Wf{c) are Wilson operators for the loops C and c, which are coaxial circles of 
radii R and r <R lying on planes separated by a length h. The correlator above defined can be 
seen as the transversal density of the colour flux generated by the source C in the fundamental 
representation and measured with the probe c. To study the flux generated by sources in a generic 
representation M, we will take Wag{C) instead of Wf{C), but keep the probe always in the funda- 
mental. 

The mean squared width of the flux tube is defined by 

The insight of [jl|] was to observe that the quantity Pf{h) can be described in the effective 
string picture by a typical Plateau problem of minimal surfaces, which can then be evaluated in the 
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infrared limit by a saddle-point approximation. More specifically, one can write 

( Wf{C)Wf{c) ) - { Wf{C) ){Wf{c) ) oc exp(-aA(/?,r,/i)) 



(1.3) 



where a is the string tension and A{R,r,h) denotes the area of the connected minimal surface with 
C and c as boundaries. Here it is assumed that the effective string is described by the Nambu-Goto 
action; it is reasonable to expect that in the infrared limit R—>oa with fixed r and h the result should 
not depend on this specific assumption. Note however that the size r of the probe cannot be too 
small, otherwise the string picture would not be valid. 

The surface with minimal area connecting the two circles is a surface of revolution about the 
symmetry axis (here chosen to coincide with the x axis). If we denote by y{x) the y >0 intersection 
of the surface with the {x,y) plane, the area is given by 



A = 271 J yy^f + \dx 



(1.4) 





By imposing the extremal condition = one finds the general solution: 

1 



i+f=yy 

ry 

express the minimal area as: 



y(x) = — cosh(o(x — xn) 
ft) 



(1.5) 



Applying the boundary conditions R = j-cosh(Oxo, r = ■^coshftj(/i — xq), it is possible to 



A = 7i(^ +R^\ll - 1 /w2/?2 - r^^i - 1 /wV 
From the boundary conditions we can express the distance as a function of (O as 

h{(o) = [Arccosh(/?ftj) — Arccosh(rft))]/ftj 



(1.6) 



(1.7) 



which is a decreasing, hence invertible, function: we can use (O as an integration variable for 
actually computing w^. Notice that since the hyperbolic cosine is always > 1 we get a minimal 
value for (O, that is a maximal allowed value for h: 

(Omin = - =^ hmax = h{l / r) 

r 

Now we can write explicitly, for the mean squared width, 

h{(oY exp[-oA{R,r,h{(o))]\h'\d(0 



(1.8) 



/r exp[-oA{R,r,h{(o))]\h'\d(0 



(1.9) 
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This quantity approaches a logarithmic curve for large R. Indeed, this can be seen by expand- 



ing Arccosh(ajr) ~ 2Iog(ft)r) for (Or ^ 1, and inserting it into ( |l.7[ ), which yelds the approximate 
solution found in rtll] : 

1 

w~-Iog(/?/r) (1.10) 
h 

The condition (Or^l becomes \og{R/r) ^ h/r, always fulfilled in the large R limit (moreover, 
r cannot be too small at fixed R and h). In this limit, a Gaussian distribution is found for the 
transversal density, whose width grows logarithmically with 7?: 



P^'^\h)ocexp - ^ aw" = -^log(/?/r) (1.11) 



OKh 



2 



log(/?/r) 



,2 



1 
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The result found here, in full agreement with the one from the conformal field theory approach 
to the problem, can be generalised to Wilson loop (C) in any representation. The worldsheet of 
the k-string can be seen as some bound state of k 1 -string worldsheets. In this case a local minimum 
solution exists, in which one of these fundamental sheets gives rise to the catenoid with the probe 
Wf{c), while the other k — I just lie flat on the loop plane. We then have immediately that the 
resulting width does not change at all from the k = I case, since 

Pl^^ih) = Plf^ih) ■ exp[7r/?2(a, - ot-i - a)] (1.12) 

2. Z4 gauge theory and duality to the Ashkin-Teller model 

The laboratory where we study the physical properties of the 2-string is a Z4 gauge model 
in three dimensions, defined by the standard plaquette action on a cubic lattice: with ci and C2 
coupling constants, and with link field G Z4 = +1,+/, — 1, — /, 



5 = £ cmiUp)+C2{Up) 



(2.1) 



With Z4 as gauge group, there exist two k-strings in the system: the fundamental string and the 2- 
string (related to the double-fundamental representation f ^ f). Moreover, since there are no more 
representations than one in the same N-ality class, the system does not exhibit any meta-stable 
string that could spoil the results at finite R. 

It is well known that this model, as any three-dimensional abelian gauge model, admits a 
spin model as its dual (see for instance [||]) and any physical property of the gauge system can be 
translated into a corresponding property of its spin dual. From a computational point of view it 
is of course much more convenient to work directly on the spin model where powerful non-local 
cluster algorithm can be applied. 

In our case the dual is a spin model with global Z4 symmetry which can be written as a 
symmetric Ashkin-Teller (AT) model [^1, i. e. two coupled, ferromagnetic, Ising models defined by 
the two-parameter action 



SaT= j8 {OxOy + TvTy) + a{OxOyZx'Cy 

<xy> '- 



(2.2) 



where Ox and are the Ising variables (a, T = ±1) associated to the site jc and the sum is over all the 
links <xy> of the dual cubic lattice. The phase diagram of this model has been studied long ago [^, 
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The global Z4 symmetry of the action is generated by the transformation a — — T ; T ^ a; the 
model has also an independent Z2 symmetry generated by a ^ T, related to the charge conjugation 
of the corresponding dual model. 

Usually the Z4 symmetry is realised by means of fields with values among the fourth roots 
of the identity = ±/,±l): the double Z2 formulation corresponds to the standard one via the 
change of fields = ^ (<7v + iTx), which preserves the above symmetries as *F ^ and the 
multiplication by /. We think it is interesting (and more purposeful for our work) to prefer the 
formulation with two Z2 coupled fields: indeed, using character expansion for the Boltzmann 
weights and resumming in terms of newly introduced (dual) link variables, it can be shown that the 
AT model is dually equivalent to two coupled Z2 gauge models with couplings a, b (Up and Vp are 
the two Z2-hke plaquette operators) 

{Ut=±l,Vi'=±l} P 

The duality is implemented by the transformation & : {cc,^) —>■ {a,b), which satisfies = 1: 

{a,b) = 



(coth j8 + tanh j3 tanh a ) (coth j3 + tanh j3 coth a ) ] 1 ^ r 1 + tanh^ j3 tanh a 

4 L tanh^jS +tanha 



2 + tanh a + coth a 



(2.4) 

2.1 Dual of a Wilson loop 

The duality transformation maps any physical observable of the gauge theory into a corre- 
sponding observable of the spin model. In particular it is well known that the Wilson loops are 
related to suitable twists of the couplings of the spin model. More specifically, let us consider a Zf^ 
spin model in 3D. Let ^ = e'^Tt/w jj^g generator of this group. A ^-twist of the link <xy> in the 
spin action is defined by the substitution j8 ^ i^*^j3 only in the selected link. Denoting with Z^xy>,k 
the spin partition function modified in this way, one can easily prove the identity 

/ u^'') ) - ^<-^>'k (2 5) 

\ '-Jp I gauge — 2 K'^--}) 
(k) 

where P is the plaquette dual to <xy> and Up is the plaquette variable in the irreducible repre- 
sentation of ZjM characterised by the integer A: = 1,2, . . . ,A'^. As a simple check, by twisting the six 
links dual to an elementary cube, one obtains the same partition function due to the invariance of 
the measure; on the gauge side, this is equivalent to measuring the product of the six plaquettes of 
the cube, which is always 1. 

In the AT model, the twist on a link can be performed regarding both Z2 components or, 
say, the a field. What results from the above recipe is that twisting only one component corre- 
sponds to operators in the fundamental representation, while twisting both is involved in the double- 
fundamental. Thus, in the AT action, we have (specifically for the frustrated link dual to the plaque- 
tte of interest) [j3 {OxGy + T^Ty) + a{axGyTxTy)] [p{-axOy + txTy) - a{axGyTx'^y)] for the funda- 
mental representation, and [jS (cj^ay + TvTy) + a(a^(7yTxTy)] [— p {OxOy + ZxTy) + a ( ay T^Ty)] for 
the double-fundamental. Combining together a suitable set of plaquettes we may build up any Wil- 
son loop or Polyakov-Polyakov correlator with k = I or k = 2 and their map into the spin model. 
Also, these are changes in the action that can be embedded in the simulation algorithm. 
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The actual measure of a plaquette operator Vp translates in magnetic terms as 

( ) gauge = { ^>-2(/3+«t.r,0a,a, (2.6) 

{Vt'^)gauge = (2.7) 

In practice, to measure a plaquette in the fundamental representation, it is a good approximation, 
at small a, use the formula (where e accounts for possible frustrations on the a component of the 
hnk): 

( ) gauge °= ( so^Oy )at + Constant (2.8) 

3. Monte Carlo analysis 

We performed a Monte Carlo analysis on the AT model, for both the fundamental and the 
double-fundamental strings, at the (confining) couphng (a,j8) = (0.0070,0.1975) (for which we 
have measured the string tensions a^a = 0.01560(1) and 0^02 = 0.0210(6)). 

On the lattice, the quantity (h) is measured by placing a square Rhy R loop W{R) on a 
plane in the desired representation and taking the probe as a plaquette operator, parallel to the loop 
and lying on its axis at a distance h. The actual measurement takes into average also the four planar 
neighbours of the plaquette in the central position, in order to enhance the signal; this operation 
does not spoil the results since we dealt with large values of 7?. 

The large loop is embedded into the update procedure as a set of frustrated links; the algorithm 
has only to measure the average value of the probe plaquette in the fundamental representation, 
since Pl^\h) oc ( )'^-^ + constant. 

The algorithm used for the analysis uses a cluster update method basically similar to the 
standard Fortuin-Kasteleyn cluster technique: each update step is composed by an update of the 
a variables using the current values of the T as a background (thus locally changing the coupling 
from j8 to j3 ± a according to the value of TvTj. on the link <xy>), followed by an update of 
the t's using the a values as background. We reached a performance of about 0.35 seconds per 
update/measurement step (as on a single Intel® Xeon 3.2 GHz 64-bit processor). 

We used a cubic lattice with side L = 80 (we found there are no finite size deviations up 
to /? ~ L/2) and measured A'^ times the plaquette operator as discussed. Configuration results have 
been then packed in groups of 25 in a binning fashion to estimate variances. We performed the 
measurements on loop sides /? = 11, 13, . . . ,41 with a statistics of 220875 measures (independently 
for each R) for the fundamental representation and 470275 for the k = 2. 

3.1 Results 

( R) 

We found that the measured values of (h) do not fit too well to a normal distribution as 
approximately expected and lead to bad width estimates, so we used the numerically integrated 
quantities instead, as in Eq. ([T^). We performed the integration by carefully choosing a cutoff 
value h,„ax- since the (indipendently measured) background value must be subtracted from the 
transverse flux density functions, the results are very sensitive to the choice of a cutoff. We took 

^max — 19. 
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Flux tube width, F vs. FF 
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By fitting tlie functions w|^j 2(^)> for R > R,nin with an appropriate distance cutting off 
non-IR contributions, to the functional form 

w^iR) = ^logR + c (3.1) 

we found that the fundamental string width, as well as the 2-string, show the expected logarithmic 
growth with the appropriate universal multiplicative factor (reduced were, respectively, 1.22 and 
2.68), but the value of c differs measurably in the two cases: ci = 0.4002(20), C2 = 0.4512(11); 
this discrepancy is probably due to some interaction between the fundamental worldsheets in the 
k = 2 case. 

Our work reinforces the numerical evidences of the broadening of the flux tube width 
extending them with high precision to the case of Z4. The main result is that this is valid also for 
non-fundamental strings: the assumption of free massless string infrared limit appears confirmed. 
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